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Abstract 

Magnetic and thermodynamical properties of itinerant-electron (metallic) 
ferromagnets described by the Hubbard model have been discussed with the 
use of the generalized Fermi-Dirac (GFD) distribution for nonextensive quan- 
tum systems. We have derived the GFD distribution within the superstatis- 
tics, which is equivalent to that obtained by the maximum-entropy method 
to the Tsallis entropy with the factorization approimation. By using the 
Hartree-Fock approximation to the electron-electron interaction in the Hub- 
bard model, we have calculated magnetic moment, energy, specific heat and 
Curie- Weiss-type spin susceptibility, as functions of the temperature and en- 
tropic index q expressing the degree of the nonextensivity: q = 1.0 corresponds 
to the Boltzmann-Gibbs statistics. It has been shown that with increasing the 
nonextensivity of \q — 1|, the temperature dependence of magnetic moment 
becomes more significant and the low-temperature electronic specific heat is 
much increased. This is attributed to enlarged Stoner excitations in the GFD 
distribution, which is elucidated by an analysis with the use of the generalized 
Sommerfeld expansion. We discuss the difference and similarity between the 
effects of the nonextensivity on metallic and insulating ferromagnets. 



PACS No.: 05.30.-d, 71.28.+d, 75.10.Lp, 75.40.Cx 



Keywords: nonextensive quantum statistics, itinerant-electron model, super- 
statistics, Hubbard model 



-'^ hideohasegawa@goo.jp 



1 Introduction 



Since Tsallis proposed the nonextensive statistics in 1988 [T] , considerable works on 
the related study have been made (for a recent review, see [2]). It is based on the 
generalized entropy (called the Tsallis entropy) which is a one-parameter generaliza- 
tion of the Boltzmann-Gibbs entropy with the entropic index q: the Tsallis entropy 
in the limit of g = 1.0 reduces to the Boltzmann-Gibbs entropy. The nonextensive 
statistics has been successfully applied to a wide class of subjects including physics, 
chemistry, information science, biology and economics [3]. Despite elegant formalism 
of the Tsallis nonextensive statistics, there are four possible methods in an evalua- 
tion of expectation values with the maximum-entropy method (MEM): (i) original 
method [1], (ii) unnormalized method [1], (iii) normalized method [5] and (iv) the 
optimal Lagrange multiplier method [6], although the four methods are equivalent 
in the sense that distributions derived in them are easily transformed each other [7] . 
A comparison among the four methods is made in Ref. |2j. 

An alternative approach to nonextensive systems is the superstatistics [H [9]. 
Complex nonextensive systems are expected to temporary and spatially fluctuate. 
It is assumed that locally the equilibrium state is described by the Boltzmann-Gibbs 
statistics, and that their global properties may be expressed by a superposition of 
them over some intensive physical quantity, e.g. the inverse temperature [8j-[ll]. 
Many applications of the concept of the superstatistics have been pointed out (for a 
recent review, see [TI]). It is, however, not clear how to obtain the mixing probability 
distribution of fluctuating parameter from first principles. This problem is currently 
controversial and some attempts in this direction have been proposed |12j-[16j. 

The nonextensive statistics has been applied not only to classical systems but 
also to quantum ones [17]-|12]. For fermion systems, the generalized Fermi-Dirac 
(GFD) distribution was derived by the asymptotic approach for \q — ll/ksT 
[T7] and by the MEM (ii) with the factorization approximation [1^. With the 
use of the MEM (ii) and MEM (iii), Refs. [HI [22] have derived the formally exact 
expression for the grand canonical partition function of nonextensive systems, which 
is expressed as a contour integral in terms of the Boltzmann-Gibbs counterpart. 
Although the exact formulation is very valuable, the actual calculation of the contour 
integral is difficult and it may be performed only in the limited cases at the moment 
|26j . Quite recently, the nonextensive quantum extension has been proposed by 
using the MEM (iv) [23]. Among the three approaches of the asymptotic [17] . 
factorization [TS] and exact methods [2TI in nonextensive quantum statistics, 
the factorization approach is the easiest way for handling physical systems [201 121] • 
The nonextensive quantum statistics has been applied to various subjects including 
black-body radiation [HI [30], Bose- Einstein condensation [Ml [251 [311 [32] , nietalhc 
[33] and superconducting materials [311 ES], spin systems [551 EZl [3S1 [3S1 HO] and 
nano-magnetism [HI |42j . 

Now we pay our attention to magnetic systems. Although there are many mag- 
netic materials, they are classified into two categories: insulating and metallic mag- 
nets. The latter are often referred to also as itinerant-electron, collective-electron or 
band magnets {metallic and itinerant- electron are interchangeably used hereafter). 
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In insulating magnets, for example, of rare-earth elements such as Gd and La, f 
electrons form the localized spin at each atomic site which yields integer magnetic 
moment in units of /is (Bohr magneton). On the contrary, in itinerant-electron 
magnets of transition metals such as Fe, Co and Ni, d electrons not only form mag- 
netic moment at each atomic site but also itinerate in crystals, by which materials 
become metallic. The magnetic moment of itinerant-electron magnets is not integer 
in units oi hb- A modern theory of magnetism has a long history over the last 
half century since the advent of quantum mechanics. Insulating magnets are well 
described by the Heisenberg model. With the Weiss molecular-field theory and more 
advanced theories for the Heisenberg model, our understanding of magnetic prop- 
erties such as magnetic structures and phase transition has been much deepen. On 
the other hand, a study of itinerant-electron magnets was initiated by Stoner [43j 
and Slater [H]. Later Hubbard proposed the so-called Hubbard model [15], which 
has been widely adopted for a study of itinerant-electron magnetism. Studies with 
the Hartree-Fock (mean-field) approximation to the electron-electron interaction in 
the Hubbard model account for the non-integral magnetic moment and the large T- 
linear coefficient of the specific heat at low temperatures, which are experimentally 
observed. 

Nonextensive statistics has been applied to insulating ferromagnets, by using the 
MEM (ii) [36][37j and the MEM (iii) [38]-|lT]. Peculiar magnetic properties observed 
in manganites are reported to be well accounted for by the nonextensive statistics 
[36]- [IQ]. The purpose of the present paper is to apply the nonextensive quantum 
statistics to itinerant-electron ferromagnets with the use of the GFD distribution 
derived within the superstatistics. The resultant GFD distribution is equivalent to 
that obtained by the MEM (ii) with the factorization approximation [TS], which is 
valid for dilute fermion gas [T^l [201 IM] • The factorized GFD distribution [TS] has 
been applied to various quantum subjects [251 EH [331 (Ml [35] because it is a 
good, practical approximation [201 |2l]. We have calculated the magnetic moment, 
energy, specific heat and Curie- Weiss-type susceptibility of itinerant-electron fer- 
romagnets described by the Hubbard model with the Hartree-Fock approximation. 
Such a calculation is worthwhile, clarifying the effect of the nonextensivity on metal- 
lic ferromagnets which is different from that on insulating counterparts. Our study 
is the first application of the nonextensive quantum method to itinerant-electron 
ferromagnets, as far as we are aware of. 

The paper is organized as follows. In Section 2, we discuss the adopted Hubbard 
model and the GFD distribution derived by the superstatistics. Analytical expres- 
sions for magnetic moments, energy, specific heat and susceptibility are presented 
with some model calculations. In Section 3, qualitative discussions on magnetic and 
thermodynamical properties are made with the use of the generalized Sommerfeld 
low-temperature expansion for physical quantities. Relevance of our calculation to 
heterogeneous magnets such as metallic spin glass and amorphous metals is dis- 
cussed. Section 4 is devoted to our conclusion. 
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2 Formulation 



2.1 Adopted model 

We have considered itinerant-electron ferromagnets described by the Hubbard model 
given by [35] 

H = ^^eon,^ + ^Y^ tij a^aj^ + ^ ^'T^U ~ ^Z^^^T - 

(T i a i,j i i 

where rijo- = aJ^Oio-, Oio- {a]^) denotes an annihilation (creation) operator of a a- 
spin electron (a =t, |) at the lattice site i, eo the intrinsic energy of atom, t^j the 
electron hopping, U the intra-atomic electron-electron interaction and B an applied 
magnetic field. We have adopted the Hartree-Fock approximation to the electron- 
electron interaction of the third term in Eq. ([T]), as given by 

Urii^riii ~ U{n^) nii + U{ni) -U{n^) (ni), (2) 

where (rzo-) denotes the average of number of electrons with spin a to be evaluated 
shortly [see Eqs. ( 1271) and ( l28l) ]. With the Hartree-Fock approximation, Eq. ([T]) 
becomes the effective one-electron Hamiltonian given by 

H ^ H^ + H^-U{n^){n^), (3) 

where 

Ha = "Yeania + ^^tijal^aja (cr =t and |), (4) 

with 

eta = eo + UriiT fJ'BB, (5) 
the minus (plus) sign in Eq. ([5]) being applied to |-spin (|-spin) electrons. 

2.2 GFD distribution within the superstatistics 

We have considered the nonextensive fermion system consisting of many clusters, 
each of which includes N particles with total energy E. Non-equilibrium or quasi- 
equilibrium states of the system are expected to temporary and spatially fluctuate. 
However, locally the equilibrium state of a given cluster is assumed to be described 
by the Boltzmann-Gibbs statistics in the superstatistics [8]-[l0]. The probability 
distribution for {n^} in a given cluster with the local temperature T is given by 

PBG0,{nk}) = ^n^"'^^''"''^"'' (6) 

where the grand-partition function H(/5) is given by 

m = Uil + e-^^^^-^^ (7) 

k 
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Here f3 = l/Zc^T, ks is the Boltzmann constant, and Uk and ek denote the number 
of states and energy, respectively, of the state k. 

After the concept of the superstatistics [8]- [10], we have assumed that the inverse 
of the temperature (/?) fluctuates and its distribution is given by the x^-distribution 
with rank n given by 



where T{x) denotes the gamma function. The average and variance of (3 are given 
by 0)^ = [3 and ((/J^)^ - (3^)/ (3^ = 2/n, respectively. The distribution averaged 
over the system with the temperature T (= l/ksP) is assumed to be given by 

POO _ _ _ 

P(AK}) = / PBGiPAnk})90)dp. (9) 

When we adopt the type-A superstatistics in which the /3 dependence of is 
neglected [11], Eqs. and yield 

p{P,{nk}) = ^ JJexpJ-^(efc (10) 
^1 k 

with 

= ll{l + exp^[-P{ek-fi)]}, (11) 

k 

q = - + (12) 
n 

where expq(x) stands for the g-exponential function defined by 

expg(x) = [l + (l-g)x]^ for 1 + (1 - g)x > 0, (13) 
= for 1 + (1 - q)x < 0. (14) 

Equation (HM expresses the cut-off properties of the g-exponential function. 
The probability of occupation of the state k by Uk particles is given by 

[l-(l-g)/?(e-/x)nfc]T^ 

+ (15) 



P{nk) 



(16) 

l + [l-(l-g)/3(e,-/i)] — 

where the factorization approximation is employed [18j. The probability of occupa- 
tion of the quantum state with the energy e^, fq{ek), is given by 

fqi^k) = ^P(nfc)nfc, (17) 

1 . . 



1 + (expj-/3(efc - /i)]) 
5 



-1 • 



Averages of the number of electrons and energy per cluster are expressed in terms 
of /g(efc) by 

N = $^/.(efc), (19) 

k 

E = 5^/,(efc)e,. (20) 

k 

The GFD distribution given in Eq. (fT8|) is equivalent to that obtained by the 
MEM (ii) with the factorization approximation for < g < 2 [18]. In the present 
paper, we will adopt the GFD distribution given by Eq. (fT8|) not only for g > 1.0 
but also q < 1.0, although it is valid for q > 1.0 within the superstatistics [Eq. (fT2|) ]. 
One of the advantages of the superstatistics is that the entropic index q is expressed 
in terms of model parameters as given by Eq. (I12p (related discussion being given 
in Section 3). 

In the limit of — > cxd where g{P) = 5{$ — 13) [Eq. ([8])] and q = 1.0 [Eq. 
f|T2l) ]. Eq. f|T8l) reduces to the conventional Fermi-Dirac distribution obtained in the 
Boltzmann-Gibbs statistics, 

/i(efe) = fBciek) = Y^r^i^- (21) 

In the zero-temperature limit of /5 — oo, both Eqs. (fT8|) and (!2T]) reduce to 

f,{ek) = /i(efc) = e(/i-efe) for /3 ^ oo, (22) 

where G(x) denotes the Heaviside function: G(x) = 1 for x > and zero otherwise. 
Equation (122|) shows that the quantum state at T = is not modified by the 
nonextensivity [21]. The derivative of /q(e) with respect to e is given by 



de {i+(expj-/?(e-/x)])-'} 



(23) 



for q 1.0. (24) 

In the high-temperature limit of (3 —>■ 0, Eq. f|T8|) becomes 

/,(e,) ~ exp^(-/5efc) for /5 ^ 0, (25) 

with /i = 0. 

Figs. [T](a) and (b) show the e dependence of /g(e) and —dfq{e)/de, respectively, 
for various q values with (3 = 1.0. When q is increased (decreased) from q = 1.0, 
the distribution at e < /i and e > is increased (decreased). The effect of q on 
the GFD distribution is more clearly seen in its derivative of dfq{e)/de, which has 
a power-law tail at |e — /i| ^ 1 for q ^ 1.0. Because of the cut-off properties of the 
g-exponential function given by Eq. (fT^ . we obtain 

/,(e) = 0.0, ^^ = at(e-^)> (^forg<l, 

/,(e) = 1.0, ^^ = at(e-/i)<-(^forg>l, (26) 
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which is clearly realized in Fig. [H 



2.3 Magnetic moment 

By using Eqs. ([3]), f|T8l) and f|T9l) with eo + Un/2 = 0, we obtain the self-consistent 
equations for the magnetic moment (m) and the number of electrons [n) per lattice 
site, given by 

m = n^-ni = J [p|(e) - Pi(e)]/g(e) de, (27) 
n = n^ + n^= [[p^{e)+p^{e)]f,{e)de, (28) 



with 

PTa(e) = Po(e± 



Um 



(29) 



where po(e) denotes the density of states and Na the number of lattice sites: the 
plus and minus signs in Eq. (129|) are applied to |- and |-spin electrons, respectively. 
From Eqs. (p71) - (!30l) . m and p are self-consistently determined as a function of T 
for given parameters of g, n and U and density of state, Po(e)- 

Bearing Fe, a typical transition-metal ferromagnet, in mind, we have performed 
model calculations with a bell-shape density of states for a single band given by 



where W denotes a half of the total bandwidth. It has been reported that U ~ 2W ~ 
5 eV for Fe |l6l HT] . Fig. [2] shows the ground-state magnetic moment as a function 
of U for n = 1.4 electrons. We have employed U/W = 1.75 which leads to magnetic 
moment of m = 0.47 fis at ksT/W = 0.0. Adopted values of n = 1.4 electrons and 
m = 0A7 fiB roughly correspond to those of Fe which has seven d electrons and the 
ground-state magnetic moment of 2.2 /i^ {i.e., 7.0/5=1.40 electron and 2.2/5=0.44 
Hb per orbital). It is noted that the U-m relation shown in Fig. [2] is valid for 
< g < 2 because /g(e) is independent of q at kBT/W = [Eq. ([22])]. We have 
solved self-consistent equations (|27|) - (|5n]) by changing q and T with the use of the 
Newton-Raphson method, which is indispensable in our calculations, in particular 
for q < 0.4 and q > 1.6 (see the appendix). 

Figs. [3] (a) and (b) show the temperature dependence of the magnetic moment, 
m, for q < 1.0 and q > 1.0, respectively. With increasing |g — 1|, the temperature 
dependence of magnetic moments becomes more significant and the Curie temper- 
ature becomes lower. This is more clearly seen in Fig. IH where Tc is plotted as a 
function of q. The q-Tc plot is almost symmetric with respect to g = 1.0 where we 
obtain the maximum value of kBTc/W = 0.143. If we adopt W ~ 2.5 eV obtained 
by the band-structure calculation for Fe [17], the calculated Curie temperature at 
q = 1.0 is Tc ~ 3500 K, while the observed Tc of Fe is 1044 K 
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2.4 Energy and Specific heat 



We calculate the energy per lattice site given by 

E = j e[p,{e)+p,{e)]Ue)de-^{n'-m% (32) 

from which the specific heat is given by 

dE dE dE dm BE dp 

dT ^ df ^ d^lr ^ ~dJIdT' ^ ' 



with 



If = -fJ^i^-^^)[P^i^)+Pli^)]^de, (34) 

^ = -^J^iM^)-PM^d^^ (35) 

^ = -J^[P^i^)+PM^de. (36) 

Analytic expressions for dm/dT and dp/dT in Eq. (!33|) are expressed by Eqs. (]A9|) - 
(lAlOp (for details see the appendix). 

Figs. [5] (a) and (b) show the temperature dependence of the specific heat C for 
q < 1.0 and q > 1.0, respectively. With increasing \q — 1.0|, the specific heat at low 
temperatures is increased and its temperature dependence is considerably modified. 



2.5 Spin susceptibility 

The spin susceptibihty is expressed by 



dm 

^ = dB^ 



from which the paramagnetic spin susceptibility is given by 



2 /or 



with 



Xo = -/ P{e)^de. (39) 

The temperature dependence of the inverse of calculated susceptibility, 1/x, for 
q < 1.0 and q > 1.0 is shown in Figs. E](a) and (b), respectively. The Curie 
temperature Tc, which is realized at 1/x = 0, is decreased with increasing \q — 1|, 
as shown in Fig. HI 

Fig. [7] shows the Curie temperature Tc as a function of U for various q, which are 
determined by the divergence of the susceptibility. The Curie temperature vanishes 
at U/W < 1.66 independently of q because the GFD does not depend on q [Eq. 
(1221) ]. The Curie temperature is lowered with increasing Iq' — 1| (Fig. H]). 
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3 Discussion 



It is possible to qualitatively elucidate the magnetic and tliermodynamical properties 
of nonextensive itinerant-electron ferromagnets presented in the preceding section, 
with the use of the generalized Sommerfeld expression of various quantities. The 
integral I including an arbitrary function 0(e) and the GFD distribution /g(e) is 
given by 

' K^)mde, (40) 

oo 

0(e)rfe + ^ c„r>("-i)(/i), (41) 



n=l 



with 



= 3- /" — ( ^ ) dx, (42) 

n\J_^ rfa; Vl + [expg(-a;)]"V 

which is valid at low temperatures. Expansion coefficients for g = 1.0 are given by 
C2 = vrVe (=1.645), C4 = 77r7360 (=1.894), and c„ = 0.0 for odd n. 

The q dependence of c„ for n = 1 to 4 is plotted in Fig. [8|, which shows the 
followings: (1) Ci and C3 are not zero for q 7^ 1.0 [33j (though a magnitude of Ci is 
small) in contrast with Ci = C3 = in the conventional Fermi-Dirac distribution, (2) 
with increasing |g — 1|, C2, jcal and C4 are much increased, and (3) the q dependence 
of C2 and C4 are almost symmetric with respect to g = 1.0 while those of Ci and C3 
are nearly anti-symmetric. The obtained q dependence of c„ may be understood as 
follows. 

When we expand /g(e) given by Eq. (fTSl) in a series of (g — 1) [24j, we obtain 
Uf-) = + + (43) 

= AW- + (44) 
Substituting Eq. (Hil) to Eq. (14T1) and using the integral by part, we obtain / given 

by 

/ = j ct^{e)h{e) de + j [2(e - ^)0(e) + (e - lif<P\e)]h{e) de. (45) 

By using Eq. fHTl) for the second term of Eq. fHS]) . we obtain 

cM - + (46) 

C2{q) = C2(l) + 0((g-l)2), (47) 

C3(.) - + (48) 

C4(g) = C4(l) + 0((g-l)2), (49) 
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where contributions of 0(g — 1) to C2 and C4 are vanishing. Thus the q dependence 
of C2 and C4 is almost symmetric with respect to q = 1.0 whereas that of ci and C3 
is nearly anti-symmetric, as Fig. [9] shows. 

Setting 0(e) = p^{e) - p|(e), 0(e) = pt(e) + Pi(e), and 0(e) = e [pt(e) + Pi(e)] in 
Eq. (14 ip with Ci = 0, we obtain (hereafter we adopt the reduced units in which 
W = ks = Hb = I) 

m(T) = m(0) + C2[p;-p'jT2 + .., (50) 
n(T) = n{Q) + C2[p\+p\]T^ + -, (51) 

^^m(T)2 + .., (52) 



E{T) = E{0) + C2[p^+p^+p{p\+p\)]T'+(^ 



4 



where p^ = Po-(p) and p^ = dp{p)/de. Simple calculations using Eqs. (j50l) -(j52l) lead 
to 

m(T) = m(0)-aT^ (53) 

C(T) = 7T, (54) 

with 

tt = C2(p'^-p'^), (55) 

7 = 2c2(pi+PT)-af/m(0). (56) 

The T^-decrease in m(T) arises from the Stoner excitations. When we take into 
account spin-wave excitations, which are neglected in the Hartree-Fock approxima- 
tion, magnetization decreases following the T'^/^ power at low temperatures. A rapid 
decrease in m{T) and a large specific heat with increasing ji? — 1| shown in Figs. [3] 
and are attributed to an enlarged C2 shown in Fig. [HI 
Setting 0(e) = dp{e)/de in Eq. fHTl) . we obtain 



Xo 



p + C2P^^'>T^ + C3p^^^T^ + -, (57) 



where p^^-* = p^^\p) {i = 2,3). The Curie temperature Tc is implicitly given by 
Uxo(Tc) — 1=0, which yields 



-C2t/p(2) 

A significant decrease in Tc with increasing Iq' — 1| shown in Fig. [His again due to 
an enlarged C2. 

For a calculation of the susceptibility at T > Tc, contributions from higher terms 
than are necessary. From Eqs. 0381) and (I37|) . the inverse of the susceptibility is 
given by 

1 ^ ^[c2p(^)(T^-Tg) + C3p(^)(T3-Tg) + -] 

X 2[p + C2p(2)T2 + C3p(3)T3 + ..] ' ^ ' 
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We note in Eq. (!59|) that if C3 = 0, the temperature dependence of 1/x becomes 
symmetric with respect to q = 1.0 because of a symmetry of C2. It is not the case 
because of the significant contribution from C3, as shown in Fig. [61 

In the present study, we have employed the GFD distribution of /g(e) given by 
Eq. ( TTSi) obtained within the superstatistics [S]-[Tn]. It is worthwhile to point out 
that the resultant GFD distribution depends on a way how the average is performed 
over the fluctuating field in the superstatistics. Indeed, if taking the average of 
fBci^k) given by Eq. fl2Tl) over the x^-distribution, we obtain 



00 



f,{ek) = I fBG{ek)giP)d(3, (60) 

g{P) dp. (61) 



00 



1 + e^(''=-'') 

The e dependence of fg{ek) calculated by numerical methods is shown in Fig. [91 We 
note that fq{ek) is rather different from fq{ek) given by Eq. (ITS]) except for q = 1.0. 
Fig. [HI clearly shows that the average over the fluctuating field in Eq. (1611) leads 
to result different from /g(e) given by Eq. (ITS]) . The chain curve in Fig. [HI will be 
explained below. 

It is noted that some applications of the nonextensive quantum statistics have 
employed the GFD distribution given by [33l [34] 

= l + (exp,|-3(«-A)|)-'' 

in place of fq{ek) in Eq. ( ITSi) . A power index q in the denominator of Eq. ( I62l) arises 
from the g-average of {0)q = Tr p'^ O in the MEM (ii), where Tr denotes the trace, 
p the density matrix and O a given operator. Such an averaging does not appear in 
either classical or quantum superstatistics [8]- [10]. The chain curve in Fig. [91 shows 
the e dependence of /g(e) for q = 1.2, which is similar to /g(e) shown by the solid 
curve. 

Magnetizations of manganites like Lao.eoYo.orCao.ssMnOs show a peculiar tem- 
perature dependence, which has been shown to arise from the cut-off properties 
in the generalized Brillouin function for nonextensive localized spin systems [36]- 
[10]. In metallic ferromagnets under consideration, the cut-off properties appear in 
the GFD distribution [Eq. (IT4])]. Because m, E and C are integrated quantities 
over the GFD distribution, effects of the cut-off properties are hardly realized in 
their temperature dependence. In both metallic and insulating ferromagnets, the 
temperature dependence of m becomes more significant with further increasing the 
nonextensivity. 

When employing Eq. (I56]) for paramagnetic metals (pj = p| = p), we obtain 
the T-linear coefficient of the specific heat given by 7 = 4c2p, which leads to the 
enhancement of the linear coefficient of the specific heat by the nonextensivity: 

7(g) C2(g) 

7 1 C2 1 
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This ratio is increased with increasing \q — 1|: it is 4.28 and 10.46 for |g — 1| = 0.5 
and 0.9, respectively [Fig. [8]. Within the superstatistics, this phenomenon may be 
interpreted as due to the effect of fluctuating [3 (for g > 1). Similar enhancements 
in the specific heat are realized in effects of spin fluctuations at low temperatures 
[ini EO] and of critical fluctuations near the transition temperatures. It has been 
reported that in some heterogeneous magnetic metals such as metallic spin glasses 
and metallic amorphous ferromagnets, contributions from electronic specific heat are 
abnormally large compared to that in normal metals. The temperature dependence 
of the specific heat and susceptibility in nano-magnets has been shown to consider- 
ably depend on the nonextensivity [SI Il2] . It would be interesting to analyze these 
materials from the view point of the nonextensive statistics. 

4 Conclusion 

By using the factorization approach to the GFD distribution [IS], we have discussed 
magnetic and thermodynamical properties of nonextensive itinerant-electron ferro- 
magnets described by the Hubbard model with the Hartree-Fock approximation. 
Our calculation has shown that with increasing the nonextensivity of |g — 1|, Stoner 
excitations is much increased, which induces a significant decrease in the magne- 
tization and a considerable increase in the specific heat at low temperatures. The 
adopted Hartree-Fock approximation has some disadvantages: it cannot well explain 
the T^/^-power law of magnetization at low temperatures, the large specific-heat 
anomaly around Tc, and the Curie- Weiss susceptibility. Nevertheless, the Hartree- 
Fock approximation has an advantage that it provides a reasonable overall descrip- 
tion for magnetic and thermodynamical properties. For a more accurate description 
of nonextensive quantum systems, it is necessary to go beyond the factorization ap- 
proximation to the GFD distribution [51] . Our study may be generalized to various 
directions: extensions to ferromagnetic (and antiferromagnetic) metals and alloys 
with more complicated structures, and calculations of various physical quantities 
such as spin waves and conductivity. We may extend our theory to include effects of 
spin fluctuations in nonextensive itinerant-electron ferromagnets with the use of the 
functional- integral method which is useful for (extensive) bulk ferromagnets [IS] [S2] • 
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Appendix: Calculations of dm/dT and dfi/dT 

The terms of dm/dT and dfi/dT in Eq. (!33|) may be derived as follows. From 
Eqs. ([27D and ([28D, we obtain 
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dT \dT 



with 



«21 ( ^ ) + 022 ( ^ ) = &2, (A2) 



(A3) 



«ii = !+(!)/ h(^) + Pi(^)]^^e, 



dm 



au = I [Me)-Pi{e)]-^de, (A4) 



«2i = l^^jy h(e)-Pi(e)]-^^e, (A5) 

«22 = J [Me) + Pi{e)]^de, (A6) 

= -^Ji^-p)[M^)-pM^de, (A7) 

&2 = -U {e-p)[p,{e)+Pi{e)]^de, (A8) 



By solving Eqs ( lAip and ( 1A2I) . we obtain 

dm (a22&i - a,i2&2) 



dT A 

dp _ {-a2ibi + 01162) 
dT ~ A 



(A9) 
(AlO) 



where A = 011022 — 012021- 

Coefficients given by Eq. ( ]A3P -( lA6l) are used also for solving Eqs. ( I271) -( l28ll by 

the Newton-Raphson method. The £th iterative solution of mi and pe of Eqs. ( 1271) 
and (128|) are given by 

(022C1 - O12C2) /A11^ 

mi = mi-i H — , (All) 

(-021C1 + 011C2) , . 

pe = pi-i H , (A12) 

with 

ci = j [p^{t) - pi{e)]fq{t) dt-mi-i, (A13) 

C2 = |h(e)+Pi(e)]/,(e)de-n, (A14) 

where the ffist terms in Eqs. ( 1A13[) and ( ]A14[) are expressed in terms of the (£ — l)th 
solutions of m£_i and pi-i. 
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Figure 1: (Color online) The energy dependence of (a) the generalized Fermi-Dirac 
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Figure 2: The U dependence of the magnetic moment m for ksT/W = 0.0, 
fisB/W = 0.0 and n = 1.4 electrons, the arrow denoting the U/W value (=1.75) 
adopted in model calculations. 
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Figure 3: (Color online) The temperature dependence of the magnetic moment m 
for (a) q < 1.0 and (b) q > 1.0 with i^bB/W = 0.0. 
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Figure 4: The Curie temperature Tc as a function of q. 



Figure 5: (Color online) The temperature dependence of the specific heat C for (a) 
q < 1.0 and (b) q > 1.0 with iibB/W = 0.0. 



Figure 6: (Color online) The temperature dependence of the inversed susceptibility 
1/x for (a) q < 1.0 and (b) q > 1.0, insets showing the enlarged plots for 0.0 < T < 
0.4. 



Figure 7: (Color online) The U dependence of the Curie temperature Tc for g = 1.0 
(the dashed curve), q = 0.2 (the chain curve) and g = 1.8 (the solid curve). 



Figure 8: (Color online) The q dependence of the coefficients c„ in the generalized 
Sommerfeld expansion [Eq. fHTj) ]: Ci (the dashed curve), C2 (the solid curve), C3 (the 
chain curve) and C4 (the double-chain curve). 



Figure 9: (Color online) The energy dependence of the GFD distributions of fq{e) 
[Eq. ([I8l)](the solid curve), fg{e) [Eq. (^] (the dashed curve) and fg{e) [Eq. (jfel)] 
(the chain curve) for q = 1.2: the three GFD distributions agree for g = 1.0 (the 
dotted curve) (see text). 
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